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Abstract. We discuss the question whether the ma tter in dusty tori around active galactic nuclei has a smooth 
or a clumpy structure. Nenkova, Ivezic & Elitzur i2002T) have argued that the lack of emission feature in the 
SEDs of type 1 AGN galaxies combined with a clear absorption feature in type 2 AGN can be explained if the 
circumnuclear dust is distributed in discrete clumps. Our aim is to verify this. We use multi-dimensional radiative 
transfer models of smooth and clumpy tori, and compare the SEDs of equivalent smooth and clumpy models. We 
find that the 10 /im emission feature of the clumpy models, when seen almost face-on, is not appreciably reduced 
compared to the equivalent smooth models. Some of the clumpy models have a weak or even absent 10 /im feature, 
but so do some of the smooth models. On the whole the SEDs of clumpy and smooth tori are similar, but some 
details are different. The absorption feature seen at edge-on inclinations appears to be less deep in the clumpy 
models than in the smooth models, and the average flux in the near-infrared regime is stronger in the clumpy 
models. Moreover, at these inclinations the clumpy models have a slightly wider SED. Whether these differences 
are unique enough to be used as a diagnostic for dumpiness of AGN tori is not yet clear. 
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1. Introduction 

According to the unification principle of active galactic nu- 
clei (AGN), the intrinsic difference between a Seyfert 1 and 
a Seyfert 2 galaxy (for the radio-quiet AGN) or between a 
narrow-line radio galaxy and a radio quasar (for radio-loud 
AGN) is m erely a questio n of orientation (Anton ucci fc 
Miller 1 19851 Barthel tl989t see review by Antonucci Il99^ . 
The leading hypothesis is that the central engine, an ac- 
creting super-massive black hole, is surrounded by a geo- 
metrically and optically thick torus of dust and gas with 
an equatorial visual optical depth much larger than unity. 
When viewed face-on, the source would look like a type 1 
active galaxy, e.g. a Seyfert 1 or a quasar, and when viewed 
edge-on, it would have the characteristics of a type 2 ac- 
tive galaxy, e.g. a Seyfert 2 or a narrow-line radio galaxy. 
From this unification model it follows that the dust in the 
circumnuclear torus emits strongly in the infrared, due to 
the irradiation by the central source. Indeed, the predicted 
infrared emission from multi-dimensional continuum ra- 
diative transfer models for such tori is reasonably well 
in agreement with the observed in frared radiation from 
such active nucl ei (Pi er & Krolik Il993t Efsthathiou & 
Rowan- Robinson 1995; Granat o, Dan ese & Francheschini 
Il997t van Bemmel & Dullemond l2003l henceforth vBD03). 
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Moreover, various other differences between Seyfert 1 / 
Seyfert 2 and between quasar / radio galaxy can be ex- 
plained in terms of such a torus model. For instance, 
the measured polarized (i.e. reflected) nuclear emission in 
type 2 sources proves that an active nucleus is present 
even t hough no direct emiss ion is observed (Antonucci & 
Miller Il985t Pier et al. ll994jh 

In spite of the success of the obscuring torus model, 
there are a number of unsolved problems with this sce- 
nario. The most troubling problem originates from the 
fact that, in order to have a hydrostatically supported ge- 
ometrically thick torus around a supermassive black hole, 
the temperature of the torus must be of the order of 
10 6 K or more. Dust in such a hot torus would not sur- 
vive long, yet dust signatures are observed to be present, 
such as a pronounced mid- to far-infrared thermal bump 
and 10 /im Si-0 stretching band of silicate in absorption 
in type 2 sources. Vario us sol utions have been proposed in 
the past. Pier & Krolik l|1992l henceforth PK92) suggested 
that radiation pressure wit hin th e torus may be enough 
to support it. Dopita et al. l)l998|) . on the other hand, put 
forward the scenario that the torus is a slowly rotating 
free-falling "envelope" , that circularizes at the centrifu- 
gal radius where it feeds the accretion disk around the 
black hole. Anot her scenario, first suggested by Krolik & 
Begelman l)l988|) , is that the torus in fact consist of a large 
number of optically thick clumps orbiting around the cen- 
tral engine and experiencing regular collisions with other 
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clumps. More recently it was shown that a nuclear star- 
burst could provide enough energy input into the torus via 
supernovae, that the torus can keep up its scale height and 
has a "sponge" like structure (Wada & Norman l2002|) . All 
of these scenarios have their strengths and problems, and 
the issue is stil l subject of debate (e.g. Vollmer, Beckert 
& Duschl l2004h . 

In spite of lack of detailed knowledge about the struc- 
ture of the torus, several studies have tried to describe its 
emission properties, using radiative transfer modeli ng; of 
smooth tori (PK02; Efsthathiou & Rowa n- Rob inson ll995l 
henceforth ERR95; Granato & Danese 1 1994 henceforth 
GD94; vBD03). They encountered two major problems in 
matching the torus emission models to the observations. 
First, many of the model spectral energy distributions 
(SEDs) were too narrow to fit the observed broad mid- 
and far-infrared SEDs in active galaxies. This can either 
be related to the presence of alternative infrared emission 
mechanisms, but it has also been shown that the adopted 
radius of the torus affects the width of the resulting SED 
(vBD03). 

A second, and still largely unsolved, issue is that 
the 10 /an silicate feature is often observed in absorp- 
tion in type 2 sources, but has never been observed in 
emission in either type 1 or 2 sources. Radiative trans- 
fer models of smooth tori tend to predict a clearly mea- 
surable 10 /im feature in emission for type 1 s ources 
(PK92, GD94, ERR95). However, Laor & Draine l|l993(l 
and vBD03 have shown that with larger grains dominat- 
ing the grain-size distribution the 10 /im feature is absent 
in type l's. Recently, Nenkova, Ivezic & Elitzur (2002, 
henceforth NIE02) proposed a different explanation: they 
suggest that clumpy tori - and only clumpy tori - natu- 
rally have these desired properties. Their claim is based on 
a model for a single clump irradiated by the central engine 
and by neighboring clumps. A statistical generalization of 
this single-clump model to a clumpy torus is made, and 
the SED computed. The clump optical depth is taken as 
a global parameter. They find that if their clump optical 
depth exceeds 60, and the typical distance between clumps 
increases proportionally to radius, then the behaviour of 
their clumpy model is in better agreement with the obser- 
vations than smooth torus models when it comes to the 
10 fj,rn feature. Relatively few clumps (typically ~ 5) are 
needed in the line of sight. In addition to this, they find 
that the SED of such a configuration is relatively wide, in 
accordance with observations. 

While the properties of the clumpy torus model of 
NIE02 are attractive, their model is highly approxima- 
tive. First of all, their single-clump model was computed 
using a 1-D radiative transfer code, even though the main 
source of irradiation of the clumps near the dust evap- 
oration radius is clearly one-sided and requires at least 
a 2-D axisymmetric approach. Secondly, their statistical 
approach to the generalization from one clump to an en- 
semble of clumps may be correct, but remains unproven. 

In this paper we take a first step toward a more self- 
consistent model and we will test the claim by NIE02 that 



infrared observations of active galaxies point to a clumpy 
torus. In order to do so, we model the clumpy torus as 
a whole, using a multi-dimensional Monte-Carlo radiative 
transfer program called RADMC. Since RADMC can only han- 
dle axisymmetric problems (i.e. 2-D problems in R and 
0), our "clumps" are in reality rings around the polar 
axis. While this setup does not constitute a realistic 3-D 
clumpy torus, it does have many of the characteristics of 
such a torus: clumps can cool by radiating in all directions, 
radiation can move freely between clumps and there are 
high density constrasts. We therefore believe that this is a 
good first step toward an understanding of the properties 
of clumpy tori. 

Our goal is to make a direct comparison between 
smooth models and clumpy models with the same global 
physical parameters. If dumpiness has a profound influ- 
ence on the SED of a torus, this comparison should yield 
distinct differences between the 10 /jm feature and overall 
width of smooth and clumpy torus models. The distribu- 
tion of the clumps is random, but on average the distri- 
bution of matter of the clumpy torus is the same as in the 
smooth torus. Following NIE02 we assume that all clumps 
have the same optical depth. 

2. Model setup 

We solve the problem of continuum radiative transfer 
through a dust density distribution around an active nu- 
cleus of luminosity L agn = 1 x 10 11 L Q . The spectral shape 
of the nuclear emission is taken to be that used in the 
models of GD94, but the precise spectral shape does not 
have a major effect on the results of model. The distri- 
bution of dusty matter around the nucleus is modeled on 
a computational grid based on spherical (polar) coordi- 
nates R, O and <&. Since our radiative transfer program 
RADMC can only handle axially symmetric density distri- 
butions, the model setup depends only on R and O. In 
addition to this, we assume mirror symmetry in the equa- 
torial plane located at 9 = 7r/2, i.e. we only model the 
domain < 8 < tt/2. The radial grid is logarithmically 
spaced, i.e. it has a constant AR/R. Such a grid ensures 
proper spatial resolution over a wide range of radii, so that 
a torus with large ratio of outer over inner radius can be 
modeled without resolution problems. The grid is lin- 
early spaced. In order to properly resolve the clumps we 
need a high spatial resolution of our grid. We have 356 
i?-grid points from the inner to the outer radius and we 
have 120 O-grid points from pole to equatorial plane. 

Our global torus setup is kept very simple. The density 
p s {R, O) for the smooth torus setup is a powerlaw function 
of i?, and is constant with O within a certain domain: 



Ps(R,6) 



p s0 {R/pc) p for abs(7r/2 - 6) < A 
for abs(7r/2 - 6) > A 



(1) 



The quantity A is the geometric thickness of the torus. 
This parameter determines the overall luminosity of the 
infrared emission of the torus, but it has only weak ef- 
fect on the shape of the SED. In this paper we shall take 
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it fixed at A = 7r/4, corresponding to an opening angle 
of 45 degrees. As inner- and outer radius we take, rather 
arbitrarily, 0.3 and 10 parsec respectively. The effect of 
varying inner and outer radius is described in detail in 
vBD03. Choosing these values differently will affect the 
width of the SED, but not so much the 10 pm feature. 
After optimizing the other model parameters, the inner 
and outer radius can be adjusted to tune the width of the 
SED and A to tune to luminosity, in order to match ob- 
servations. Many of the additional effects of the geometry 
of the torus have been described already by vBD03 and by 
others (e.g. ERR95; GD04), so we do not need to repeat 
all of them here. 

We present four smooth torus and sixteen clumpy 
models in this paper. For each smooth model there are 
four clumpy models with the same global physical param- 
eters. Two of the clumpy models have 40 clumps, and two 
have 20 clumps. The only difference between the pairs of 
clumpy models with identical number of clumps are the 
random positions of the clumps in i?, G. The smooth mod- 
els will be used as benchmarks against which the clumpy 
models can be compared. In this paper we will focus on 
the effect of the radial powcrlaw p of the density distribu- 
tion (p (x R p ), and study the effect of thermal decoupling 
between silicate and graphite grains. 



2.1. Smooth torus description 



2.2. Clumpy torus description 

An overview of the clumpy models is given in Table |21 
For the clumpy torus models we start from the smooth 
torus models SI- • -S4, and contract the matter into dis- 
crete annular clumps randomly positioned on the R,Q 
computational grid. The random positions are distributed 
such that on average the density is the same as that of the 
equivalent smooth torus. Following NIE02 we take the op- 
tical depth of the clumps to be a global parameter of the 
model. Another global parameter is the relative size of the 
clumps compared to R, i.e. a = size/R. This means that 
the size of the clumps scales with distance from the black 
hole. For the models presented in this paper we take this 
constant to be a — 0.025. Ideally we would wish to model 
smaller clumps (and more of them), but technical limita- 
tions of the resolution of our computational domain also 
limit the minimum size of our clumps, since it is important 
that all clumps are well resolved by the grid. 

We present two pairs of clumpy models to match each 
smooth torus model, the clumpy models are numbered 
accordingly, i.e. CI equals SI, etc. Between the pairs of 
clumpy models, only the number of clumps is varied, which 
we denote with suffices a and b, the a-series always having 
40 clumps and the b-series having 20. Within the pairs we 
vary the random distribution of the clumps, allowing us to 
also study the effect of randomness on the resulting SED. 
This is denoted with the number following suffix a or b. 

The matter within each individual clumps is dis- 
tributed as follows: 



An overview of the smooth models is given in Table ^ 
All models assume a 50% carbon and 50% graphite mix- 
ture for the dust. For the silicate opac ity we use the op- 
tical constants of Laor & Draine (^993). The optical con- 
stan ts for amorphous carbon were taken from Preibisch et 
al. |l993). All models have a total dust mass in the torus 
of 2 x 10 6 M Q . We vary p between and -1. We also vary 
the dust properties: on the one hand we use thermally 
coupled 0.25 /jm sized silicate and carbon grains without 
scattering (the scattering opacity taken to be zero); on 
the other hand we use a distribution of Galactic (MRN) 
dust between 0.005 and 1 /im with a treatment of scatter- 
ing (albeit in isotropic approximation). A standard MRN 
distribution only extends up to 0.25 /im, but we chose to 
extend it to 1 pm to maximize any possible effects of such 
a distribution compared to the single-grain-size models. It 
should be noted that we have only introduced a thermal 
decoupling between the carbon and silicate grains, but 
not between the different sizes of the MRN distribution, 
since the latter would be very computationally expensive, 
in particular for the high spatial grid resolution required 
for the clumpy models described below. We do not expect 
this to have much effect, since grains up to 1 pm size all 
have the same 10 pm feature shape. But future modeling 
will have to verify this. 



, . / (R/R.-1) 2 (e-6i) 2 \ 
Pci (R,e)=p c(H eiq>{-^-^ — jT^) > 

(2) 

where i stands for clump number i. 
2.3. SED generation 

Once the clumps are put onto the computational do- 
main, and the density distribution p(R, 0) is set up, 
the continuum radiative transfer problem is solved us- 
ing a Mo nte-C arlo program called RADMC (Dullemond & 
Dominik 2004), which uses an improve d ver sion of the 
original algorithm of Bjorkman & Wood l(200ll) . This pro- 
gram solves the transport of continuum radiation and the 
local thermal equilibrium of the dust grains, thus obtain- 
ing the temperature of the dust everywhere on the grid. 
The SED at inclinations of 20,70, and 90 degrees are then 
computed using a ray tracing code. We choose i = 20° to 
be the representative inclination for type 1 (face-on) AGN 
since with our sharp inner edge setup a perfect i = 0° incli- 
nation would produce an artificially low near-infrared flux 
because one would look perfectly along the inner edge. 
The representative inclination for type 2 (edge-on) AGN 
is chosen to be i = 90°. 
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Model No. 


P 


GSD 


TD 


scat 


Wi 


W 2 


Si 


s 2 


1 


SI 





0.25 






1.17 


0.44 


0.314 


-2.894 


0.110 


S2 


-1 


0.25 






1.13 


0.54 


-0.016 


-1.467 


0.068 


S3 





MRN 


V 


V 


1.26 


0.41 


0.137 


-2.025 


0.090 


S4 


-1 


MRN 


V 


V 


1.13 


0.51 


-0.083 


-1.906 


0.060 



Table 1. Overview of parameters and results of the smooth models. Columns from left to right: model number, 
powerlaw index p for p(R) oc R p , grain size distribution (either 0.25/im or MRN distribution), temperature decoupling 
between silicate and carbon grains, inclusion of scattering opacity, the resulting SED widths Wi (for i = 20°) and W2 
(for i — 90°), the resulting 10 /im feature strengths <Si (for i = 20°) and S2 (for i = 90°) and the resulting anisotropy 
parameter X. See text for definition of Wi , W2 , Si , S% and X. 



3. Results and analysis 

In this section we present the results of the model calcu- 
lations. We present figures of the SEDs at different incli- 
nation angles. We analyze our results using quantitative 
numbers for the width W of the SED (W% for face-on and 
W2 for edge-on), the strength S of the 10 /im feature (<Si 
for face-on and 1S2 for edge-on) and the anisotropy param- 
eter X. Following PK92 and GD94 we define the width W 
(i.e. Wi and W2) of the SED as the 10-log of the frequency 
range in which the spectrum is more than 1/3 of its peak 
value. In vF v ; for a pure blackbody spectrum this value 
is W= 0.686, but observations indicate that the SED of 
active galaxies is generally much wider. 

The feature strength S (i.e. <Si and <S 2 ) is defined as 
the e-log of the peak-over-continuum ratio of the feature 
for face-on inclinations. Here an e-log is used, allowing 
a direct comparison to previous studies (GD94, Laor & 
Draine 1993). Following GD94 the continuum is defined 
by a powerlaw connecting the fluxes at 6.8 and 13.9 /im. A 
positive value of S means the 10 /im feature is in emission, 
a negative value means absorption. 

Finally we define the isotropy parameter X as the linear 
ratio of the total integrated infrared flux at 90° inclina- 
tion over the total integrated infrared flux at 20° inclina- 
tion. This implies that for larger values of X there is more 
isotropy, X= 1 indicating perfect isotropy. 

3.1. Smooth models 

The results of the smooth torus models SI- • • S4 are shown 
in Fig.^ We show the density distribution on the left, the 
total SED at different inclination angles in the middle, and 
on the right a zoom-in on the 10 /im region for face-on 
inclination. 

For model SI we have taken a constant density 
throughout the torus, because according to GD94 such 
models fit better the observed SEDs of AGN. We find 
that the SED has a width of Wi= 1.17 (see table HJ. The 
model produces a quite strong 10 micron emission feature 
for face-on tori (S%= 0.314), which is inconsistent with 
observations. This strong emission feature is not entirely 
surprising, because the torus is not very optically thick. 
At R = lpc the 10 /*m vertical optical depth through the 
torus is 2.9, and just shortward of the feature at 7 /im 
it is 1.0. In fact, the nuclear radiation impinging on the 



inner rim of the torus only reaches its r = 1 surface (at 
A ~ 1/im) at a radius of about R = 0.42 pc which is about 
1.5 times the radius of the inner rim. Therefore, the in- 
ner regions of the torus are relatively optically thin, and 
expected to produce a strong emission feature. 

Model S2 is identical to SI except for changing p = — 1. 
Effectively, this puts more mass in the innermost regions 
of the disk, making these regions more optically thick. As 
a result of this, the 10 micron emission feature virtually 
disappears (S± = —0.016). The SED also becomes a tiny 
bit narrower (Wi = 1.13). These effects are similar to what 
was already found by PK92: compact tori with very high 
optical depth can have a very weak emission feature, but 
produce narrower SEDs (see also discussions in GD94). In 
the present case our torus is still relatively large and hence 
not really compact. The large outer radius ensures that 
there is a large reservoir of cool dust, which is relatively 
unaffected by the change in p. This is why the width is 
not strongly affected when varying p. The 10 /im feature, 
on the other hand, comes from the warm inner regions 
(i.e. small radius), which are very optically thick, hence 
the disappearance of the 10 /im feature. 

It should be noted, however, that the disappearance of 
the 10 /im feature for highly optically thick inner regions 
of the torus is not obvious. It depends strongly on the 
geometry of these inner regions. In our conical torus model 
the edges are very sharp and straight. There is only one 
surface of the torus that is directly irradiated, which is the 
inner rim. When this irradiated hot inner rim is seen under 
an inclination of only i — 20°, as in the figures, then the 
line of sight toward a surface clement of the rim has a high 
inclination (i ~ 70°) with respect to the normal vector on 
the rims surface, i.e. almost parallel to the surface. This 
weakens the emission feature. Moreover, the emission from 
the near side of the rim is partly re-absorbed again by the 
rim material itself. Additionally, a significant fraction of 
the 10 /im flux comes from larger radii than the inner 
rim radius. That emission has (if anything) an absorption 
feature. The end effect is that if our torus is optically thick 
near the inner rim, the emission feature is very weak. The 
total flux from the torus (including the regions at larger 
radii) may even have a slight absorption feature, even for 
i = 20°, as can be seen in model S2 with Si — —0.016. 

Models S3 and S4 are like model SI and S2 respec- 
tively, but now with decoupling the temperatures for 
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Model No. 


P 


■^V clump 


GSD 


TD 


scat 


'clump 


Wi 


W 2 


Si 


s 2 


1 


Cl_al 





40 


0.25 






26 


0.91 


0.57 


0.347 


-1.585 


0.176 


Cl_a2 





40 


0.25 






26 


0.91 


0.63 


0.377 


-0.486 


0.301 


Cl.bl 





20 


0.25 






53 


0.83 


0.60 


0.330 


-0.396 


0.250 


Cl_b2 





20 


0.25 






51 


0.77 


0.66 


0.433 


0.039 


0.294 


C2_al 


-1 


40 


0.25 






60 


0.96 


0.63 


0.170 


0.068 


0.089 


C2_a2 


-1 


40 


0.25 






60 


1.05 


0.63 


0.117 


-0.880 


0.098 


C2_bl 


-1 


20 


0.25 






105 


1.00 


0.66 


0.318 


-0.373 


0.126 


C2_b2 


-1 


20 


0.25 






115 


1.17 


0.92 


0.150 


-0.464 


0.157 


C3_al 





40 


MRN 


V 


V 


53 


1.00 


0.55 


0.260 


-0.953 


0.161 


C3_a2 





40 


MRN 


V 


V 


53 


0.92 


0.63 


0.263 


-0.381 


0.280 


C3_bl 





20 


MRN 




V 


109 


0.92 


0.60 


0.256 


-0.257 


0.240 


C3_b2 





20 


MRN 


V 


v 7 


104 


0.83 


0.66 


0.320 


-0.025 


0.289 


C4_al 


-1 


40 


MRN 


V 


V 


123 


1.05 


0.63 


0.131 


-0.190 


0.086 


C4_a2 


-1 


40 


MRN 


V 


V 


123 


1.05 


0.75 


0.059 


-0.703 


0.085 


C4_bl 


-1 


20 


MRN 


V 


V 


214 


1.09 


0.66 


0.245 


-0.336 


0.121 


C4_b2 


-1 


20 


MRN 


si 


V 


235 


1.17 


0.71 


0.091 


-0.360 


0.138 



Table 2. Overview of the model parameters and results of the clumpy models. Columns from left to right: model 
number, powerlaw index p for p(R) oc R p , number of clumps, grain size distribution (either 0.25/im or MRN distribu- 
tion), temperature decoupling between silicate and carbon grains, inclusion of scattering opacity, optical depth of the 
clumps, the resulting width VVi of the face-on SED, W2 for the edge-on SED, the resulting 10 /im feature strength Si 
for the face-on SED and 6>2 for the edge-on SED and the resulting anisotropy parameter X. See text for definition of 
Wi, W2, <Si, S2 and X. The optical depth of the clumps follows from the number of clumps, the distribution of clumps 
and the total mass of the torus, and is therefore not an independent parameter of the model, hence the separated 
column for the optical depth. 



graphite and silicate grains, introducing an MRN grain 
size distribution and including isotropic scattering. We 
find that the 10 fj,m emission feature is significantly weak- 
ened in S3 compared to model SI (6>i= 0.137 instead of 
<Si= 0.314), but still clearly present. The reason for the 
weakening is that the graphite has a higher opacity at 
visual wavelengths than silicate. Therefore the graphite 
becomes hotter than the silicate, and produces more con- 
tinuum emission in the 10 /im wavelength region. In model 
S4 we see the silicate feature slightly in absorption (Si = 
-0.083) for i = 20° again. The width of the SED is vir- 
tually unaffected by changing to a proper dust grain size 
distribution. If anything, S3 is somewhat wider than SI, 
but no difference is found between S4 and S2. The same 
goes for the anisotropy: S3 is slightly more anisotropic 
than SI, but S4 and S2 do not differ. 

3.2. Clumpy models 

The clumpy models CI- • • C4 are the clumpy generaliza- 
tions of models S1---S4. The results for the quantities 
Wi, W 2 , Si, S 2 and X are given in Table H The SEDs are 
shown in Fig. [21 for models CI and C2, and in Fig. [3] for 
models C3 and C4. Again, on the left the density distribu- 
tion, middle the total SED at different inclination angles, 
and on the right a zoom-in on the 10 /im region. 

Figs. 12131 show that with the still relatively large size 
of the clumps (a = 0.025) and the number of clumps used 
in these models, the distance between clumps is not al- 
ways very much larger than the size of the clumps them- 
selves. The NIE02 definition of a clumpy torus requires 



that the clumps are very small compared to the mean 
free path between the clumps. In this respect our model 
is not a true representation of clumpy tori as defined by 
NIE02. However, as mentioned above, technical limita- 
tions prevent us at present from modeling smaller and 
more clumps. We use the a- and b-series to understand 
the effect of the inter-clump distance, which is larger in 
the b-series. 

The values of Wi do not seem to follow a major 
trend, except perhaps that models C2 and C4 produce 
marginally broader SEDs than models CI and C3. This 
is the opposite effect of what we observe in the smooth 
models, where S2 and S4 generate narrower SEDs than 
SI and S3. The isotropy X has a somewhat stronger trend 
with the p = —1 ('b') models being more isotropic than 
the p = ('a') models, and the models CI and C3 being 
more isotropic than models C2 and C4. 

The random location of the clumps can have a reason- 
ably strong effect on the values of Wi , W2 and X and even 
more so on the values of Si and 1S2. The reason for these 
trends is more difficult to understand than for the smooth 
models. 

3.3. Comparison between smooth and clumpy models 

The face-on SED of models CI are clearly less wide (0.77 < 
Wi < 0.91) than for model SI (Wi= 1.17), despite the 
same global parameters. The same is true for model C3 
compared to S3. The difference in Wi between C2-S2 and 
C4-S4 is less pronounced, but still present. For edge-on 
inclinations the clumpy models have clearly wider SED 
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Fig. 1. The results of the smooth torus models. Each row represents one model. The leftmost column represents the 
density distribution. The middle column is the SED of the system at three inclinations: 20°, 70° and 90° (where 20° 
is effectively face-on, i.e. corresponding to a perfect type 1 source and 90° is edge-on, i.e. corresponding to a perfect 
type 2 source.). The vertical axis represents 10 \og(vF u ) as seen at a distance of d = 10 7 pc. The SED includes the 
emission from the torus as well as the emission from the illuminating star, though the latter is only marginally visible 
in this figure. The right column is a linear plot of the face-on (i=20°) spectrum in the 10 /im regime, again vF v at 
d = 10 7 pc. 



than the smooth models with e.g. model SI having W2= 
0.44 and the CI models having W2 ranging between 0.57 
and 0.66. 

The clumpy models have a more isotropic emission, as 
can be seen from the values of I. For instance, the model 
S3 has 1= 0.09 while the models C3 have values ranging 
between X= 0.16 and 1= 0.29. 

The strength of the 10 /im emission feature for the 
i — 20° inclination does not seem to be decreased in the 
clumpy models. In fact, on average the value of S\ is some- 
what larger for the clumpy models than for the smooth 
models. It is striking that while the S2 and S4 models 
have the feature slightly in absorption (S\= —0.016 and 
<Si= —0.083 respectively), the corresponding clumpy mod- 
els C2_kx and C4_xx have it consistently in emission: for 
the case C2_bl even strongly in emission {S\= 0.32). So 
the dumpiness apparently does little to suppress the emis- 
sion feature. It does, however, appear to weaken the ab- 



sorbtion feature for the high inclinations (type 2 AGN). 
For example, the model S3 has £2= —2.03, which is a 
strong absorption feature, while the models C3_xx have 
values ranging between —0.95 and —0.03. 

4. Discussion 

4.1. Silicate feature strength in clumpy tori 

The results of our models indicate that dumpiness does 
not seem to have the effect of suppressing the 10 /im sil- 
icate feature. In fact, the smooth models appear to do 
better in this respect than the clumpy models. Some of 
our clumpy models do have a suppressed 10 /im feature 
(e.g. model C4_a2 with Si= 0.06), but that seems to be 
strongly influenced by the random location of the clumps 
(model C4_al has Si— 0.13). One reason why the C2 and 
C4 clumpy tori have on average stronger 10 /zm emission 
features in their face-on SEDs than the smooth models 
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Fig. 2. The results of the clumpy torus models CI and C2. See Fig. for explanation. In contrast to Fig.^ however, 
the left panel shows p ■ R instead of p, so that all clumps have the same grey depth. 
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Fig. 3. As Fig. |3 but for models C3 and C4. 
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S2 and S4 (which have no emission feature) is that the 
chimps have a 'fluffy' photosphere while the smooth mod- 
els have sharp edges. Emission features are best produced 
if a photosphere is irradiated under a reasonably small 
angle, which is the case for a significant portion of the 
fluffy photosphere of the clumps. Another reason is that 
in the smooth models a part of the emission from the in- 
ner rim (in particular the emission from the near side) is 
re-absorped by the torus, producing an absorption feature 
there. For the clumpy tori this effect is weaker. 

The weaker absorption feature for edge-on clumpy tori 
compared to the edge-on smooth tori appears to be due 
to the fact that for clumpy tori one can see at least partly 
between the clumps deeper into the torus where the 10 /im 
emission is produced. Moreover, each of the clumps in the 
outer regions of the torus, having a higher density than the 
equivalent smooth torus, either blocks the light from the 
inner regions entirely (if it is in the line of sight) or does 
not block any light (if it is out of the line of sight). Only in 
few cases does one look just through the fluffy photosphere 
of the clump, causing a 10 ixm absorption feature. In other 
words: the line-of-sight extinction caused by a collection 
of clumps is more 'grey' than that of smoothly distributed 
matter. 

It should be kept in mind that our results may depend 
on the opacity we use. To allow a direct comparison with 
the NIE02 models, we have chosen a Galactic grain size 
distribution and do not use larger grains than 1 /im. Grain 
sizes and distribution can have a profound ef fect o n the 
silicate feature, as shown by Laor & Draine ll 1993^1 and 
vBD03. There is convincing observational evidence that 
the dust in active g alaxie s does not have Galactic proper- 
ties (Maiolino et al. l200l|) . Also, in radio galaxy NGC 4261 
the near-infrared colours of the observed 300 pc scalee disk 
cann ot be modeled with standard Galactic dust (Martel 
et al. l2000fl . Therefore, it is quite conceivable that the sil- 
icate feature is in reality weakened by opacity effects. On 
the other hand, it is also important to ask the question 
whether our simplification of thermally coupling the dif- 
ferent grain sizes would perhaps artificially suppress the 
feature. 

4.2. Width of the SEDs 

As mentioned above, the width of the face-on SED Wi is 
slightly smaller for the clumpy models than for the smooth 
models, in particular for the p = models (SI, S3, CI, C3). 
This effect can be explained by the fact that for p = 
the typical number of clumps per <51ogi? is not constant. 
For these models the clumps are typically more concen- 
trated to the outer regions of the torus. Since the number 
of clumps is relatively moderate, the chance is then high 
that there will be few (if any) clumps in the inner regions. 
The real inner radius of the torus is determined by the lo- 
cation of the clumps. Therefore this deficiency of clumps 
effectively shifts the inner radius of the clumpy torus out- 
ward. If models would be made with p < — 1, the same 



reasoning would be applicable for the outer radius of the 
torus. Since the width of the face-on SED depends largely 
on the ratio of the outer to the inner radius (see e.g. GD94, 
vBD03), the effect of this shift will be to narrow the SED. 
For the case of p > —1 the short- wavelength part of the 
SED will be suppressed and for the case of p < —1 the 
long-wavelength part of the SED will be suppressed. 

On the other hand, for the edge-on SED one would 
expect that clumpy tori have less absorption of the near-IR 
flux because emission from the inner regions of the torus 
can travel between clumps toward the edge-on observer. 
This effect is indeed seen in the models, which have larger 
W2 values for the clumpy models than for the smooth 
models. 



4.3. Isotropy of infrared emission 

Like a smooth torus, the bolometric infrared flux of a 
clumpy torus is stronger in the polar direction than in 
equatorial directions. This is because in both cases the 
emission from the hotter inner regions cannot be directly 
observed for edge-on systems due to the obscuration by 
the cooler outer regions. This anisotropy of the infrared 
emission is typical for all disk- or torus-like configura- 
tions. Our models of clumpy tori show that their SEDs 
are generally more isotropic than their smooth counter- 
parts. Some models show this stronger than others. In 
particular we find this effect to be very strong for models 
CI compared to model SI, while the effect is much less 
pronounced for models C2 compared to model S2. The 
reason for the increased isotropy of clumpy models is that 
radiation can freely move in between clumps. In a way, the 
intcr-clump distance acts as a new kind of mean-free-path, 
and since there are not so many clumps in the model, the 
new effective optical depth for the clumpy torus is there- 
fore smaller (even though the actual optical depth along 
individual lines-of-sight may be still high). This effect is 
also reflected in the fact that the edge-on mid-infrared flux 
for the clumpy models is significantly higher than for the 
smooth models. 

The increased isotropy in clumpy models CI and C3 
compared to C2 and C4 is for a large part due to the fact 
that there were not so many clumps at small radii, and 
therefore the inner radius was effectively moved outward 
(as discussed in subsection l4.2l above'l. Since the total mass 
of the torus was kept constant (at 10 6 M & ) and the num- 
ber of clumps was kept the same (A" =40 resp. A" =20), 
the optical depth of the clumps is lower for the models CI 
and C3 (p = 0) than for the models C2 and C4 (p = -1). 
A lower optical depth increases the isotropy of the torus 
(a perfectly optically thin torus being perfectly isotropic). 
The generally higher isotropy of the N =20 models versus 
the N =40 models (in spite of the lower clump optical 
depth of the latter) is because for fewer clumps the inter- 
clump distance is larger and the effective optical depth of 
the torus is decreased. 
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4.4. A criterion for dumpiness? 

It is interesting to ask if, instead of the 10 fim silicate emis- 
sion feature, some of the other aspects of the infrared spec- 
tra of tori could give clear indications for dumpiness. For 
example, for edge-on inclinations (type 2 AGN) the width 
of the SED, the mid- over far- infrared flux ratio and the 
depth of the absorption feature do show some trend with 
dumpiness: the clumpy models having a slightly wider 
SED, higher mid-infrared flux and shallower absorption 
feature. Unfortunately the current study is not compre- 
hensive enough to assure that these effects could not be 
reproduced by other properties of the torus, such as the 
opacity properties or the torus geometry. One has to scan 
a large parameter space of smooth models to make sure 
that some property of the spectra of clumpy tori is really 
unique to clumpy tori. This is, however, beyond the scope 
of this paper. 

4.5. The nature of clumps and dumpiness 

The models presented in this paper are meant to verify 
what is the effect of dumpiness on the SED of torus models 
for active galaxies. Yet how representative are our models 
for clumpy tori? Our clumps are not real 3-D clumps due 
to our 2-D approximation, but we have argued that this 
should not have a major effect on our conclusions. Aside 
from the 2-D issue, is our description of the clumps realis- 
tic? Very little is known about the structure of the dusty 
circumnuclear matter in active galactic nuclei. Arguments 
for dumpiness have so far been rather indirect, but it 
seems reasonable to assume that the circumnuclear mat- 
ter is distributed in an irregular and chaotic way rather 
than in a smooth and ordered manner. 

The kind of dumpiness would depend much on the 
mechanism causing the clu mps. According to the model of 
Krolik & Begelman ( 1988) these clumps are individual dy- 
namically independent objects orbiting the black hole, and 
experiencing regular semi-elastic collisions. These clumps 
must be very compact, self-gravitating, and must be sup- 
ported by strong interal magnetic fields to provide suf- 
ficient elasticity upon collisions with other clumps. On 
the other hand, a supersonically turbul ent m edium of the 
kind described by Wada Sz Norman lj2002j) would pro- 
duce filamentary (sponge) structures rather than isolated 
clumps. These different structures of dumpiness may have 
very different infrared emission properties. For instance, 
a filamctary medium is likely to have more matter in a 
marginally optically thin state than a medium consist- 
ing of very compact clumps. Since an emission feature 
comes from marginally optically thin regions, such an 
emission feature is expected to be stronger for the fil- 
ametary medium than for the compact clumpy medium. 

In the light of this, it is interesting to question what 
the effect is of the "fluffiness" of the clumps in our simula- 
tion. In contrast to NIE02 we assume our clumps to have 
a Gaussian density profile. If we would take constant den- 
sity clumps with a sharp edge, like NIE02, these clumps 



may have less marginally optically thin material at their 
surface, perhaps suppressing thereby the 10 /im emission 
feature where our models exhibit this feature clearly in 
emission. On the other hand, even for a perfectly sharp 
edge of an optically thick clump, it is not guaranteed that 
the feature vanishes because the surface of such a clump 
may be super-heated with respect to the clump interior 
by the irradiation, yielding a hot optically thin emission- 
feature-producing layer similar to what was descri bed fo r 
flared circumstellar disks by Chiang & Goldreich 
It is, unfortunately, not possible for us to investigate this 
with our current models because this would require us 
to increase the resolution of our computational grid by 
a large factor in order to make sure to sample the photo- 
sphere of the clumps properly. This would be prohibitively 
computationally expensive at present. We can therefore 
only draw conclusions about clumpy media with Gaussian 
clumps. 

Another issue related to this is the assumption, made 
by both NIE02 and ourselves, that all the clumps have 
equal optical depth. According to Krolik & Begelman 
l|l988|) the clumps in a circumnuclear torus get regularly 
tidally disrupted, forming smaller clumps which subse- 
quently merge to form bigger ones. In effect an equilib- 
rium distribution of clump sizes will result, with clumps 
of various sizes (and optical depths) coexisting within the 
same torus. Some of these clumps may easily be optically 
thin, or at least have low optical depth. Such clumps may 
again provide a reservoir of marginally optically thin ma- 
terial which could produce a 10 /im feature in emission. 
How strong this effect will be depends on the equilibrium 
distribution function of clump sizes. 

Finally it is important to mention that due to technical 
limitations we could only model a rather limited number 
of clumps, each with a rather large size. It cannot be ex- 
cluded that some results may change if we would be able 
to model problems with a much higher number of clumps, 
all of which being much smaller than we have assumed in 
the models we presented here. It is hard to estimate how 
big these effects are. In the limit of increasingly many 
ever smaller clumps, while keeping the total mass of the 
torus and the clump filling factor (the average number of 
clumps along the line of sight toward the center) constant, 
the clump optical depth eventually drops below unity. In 
this case the SED would become identical to that of the 
smooth version. If one, on the other hand, keeps the op- 
tical depth of the clumps constant while increasing the 
number of clumps, the average number of clumps along 
the line of sight drops below unity, which would be against 
the whole idea of obscuring circumnuclear tori. One would 
have to increase the total mass of the torus to compensate 
for this. Since it poses technical problems to model much 
smaller clumps than we have done in this paper, we can- 
not be certain what the effect of such an increase of the 
number of clumps would be. But by comparing the mod- 
els with N = 40 and N — 20 (the a-series to the b-series) 
we find very little differences, so we expect this to remain 
this way for very high N. 



Dullemond & van Bemmel: Clumpy tori around AGN 



11 



5. Conclusion 

We present the first global simulations of clumpy tori 
around AGN using a axisymmetric, multi-dimensional ra- 
diative transfer model. From our analysis and comparison 
between smooth and clumpy tori models we conclude that 
the 10 ^m feature can both be strengthened and weakened 
when dumpiness is introduced. The width of the SED is 
largely determined by the inner and outer radius and the 
main effect of dumpiness is to increase the effective in- 
ner radius and/or decrease the effective inner/outer ra- 
dius due to statistical fluctuations in the positioning of 
the clumps. This results in a slightly narrower face-on 
SED. The edge-on SED for clumpy tori is slightly wider 
than for smooth tori, because radiation can move freely 
between the clumps and emission from the inner regions 
of the clumpy torus can more easily reach the observer 
even for edge-on inclinations. We find that the isotropy of 
the infrared emission is significantly affected by dumpi- 
ness for similar reasons. Unfortunately this is a quantity 
that cannot be directly observed in an individual source. 
It requires studies of large samples of sources at different 
inclinations with comparable physical properties. 

We do confirm that for the particular parameters of 
dumpiness mentioned in NIE02 the 10 ^m feature can 
be rather weak, but this is even more pronounced for a 
smooth torus with the same global parameters and aver- 
age density. A stronger effect is the depth of the 10 yum 
absorption feature for edge-on tori: for clumpy tori it is 
clearly less deep than for smooth tori. 

Although we use clumps with a slightly different struc- 
ture than NIE02, and although the clumps in our 2-D 
models are annuli around the symmetry axis instead of 
real 3-D clumps, we believe that our models produce at 
least qualitatively the correct results for clumpy media. 
We therefore cast doubt on the idea that the properties 
of the 10 iim feature of type 1 and type 2 active galax- 
ies point unequivocally to clumpy tori. It should be clear, 
though, that we do not claim that the circumnuclear mat- 
ter is smooth. We merely call for caution in interpreting 
the properties of the SEDs in the context of dumpiness of 
the dusty tori in active galaxies. 
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